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Attention  is  called  to  the  fact  that  the  Wideband  satellite  data  have 
consistently  shown  phase  spectral  density  functions  that  fall  off  less  steeply 
than  f~3,  whereas  recent  barium  data  from  STRESS  show  more  steeply  sloped 
spectra.  In  light  of  the  impact  on  propagation  modeling,  it  is  important  that 
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tency  of  the  data  analysis  procedures  used  to  estimate  the  phase  power- law 
index. 


UNCLASSIFIED 


CONTENTS 


LIST  OF  ILLUSTRATIONS .  2 

I  INTRODUCTION  .  3 

II  THE  MUTUAL  COHERENCE  FUNCTION .  7 

III  THE  PHASE  STRUCTURE  FUNCTION  .  10 

IV  THE  TAYLOR  SERIES  APPROXIMATION .  19 

V  SUMMARY  AND  DISCUSSION .  22 

REFERENCES .  24 

APPENDICES 


A  THE  TAYLOR  SERIES  APPROXIMATION  TO  THE  PHASE  STRUCTURE 
FUNCTION . . 

B  THE  SMALL  q  APPROXIMATION . . 

O  . 


m. 


1  mrcssiON  for  Lrf 

NTIS 

White  Section 

ODC 

Buff  Section  □ 

UNANNOUNCED 

□ 

JUSTIFICATION - 

|  OQTRIBIITION/AVMUBUn  MB  J 

lOist  WAIL  if*/* 

/} 

.| 

1 


a  .  -  -viir.  - 


ILLUSTRATIONS 


1  Plot  of  Normalized  Phase  Structure  Function  Showing 

Effect  of  c  for  Shallowly  Sloped  Phase  SDF  ••••••••••  13 

2  Plot  of  Normalized  Phase  Structure  Function  Showing 

Effect  of  c  for  Steeply  Sloped  Phase  SDF  •  •  •  . .  14 

3  Plot  of  Normalized  Structure  Function  Showing 

Dependence  on  Spectral  Slope  .  15 

4  Plot  of  Small  q0  Approximation  to  Phase  Structure 
Function  (valid  for  v  <  1*5)  Superimposed  on  Exact 

Curves . . . .  •  16 

5  Plot  of  Quadratic  Approximation  to  Phase  Structure 

Function  Superimposed  on  Exact  Curves  .  21 


i 


2 


I  Li*.} 


I  INTRODUCTION 


In  this  report  the  ramifications  of  some  new  theoretical  results 
for  modeling  late-time,  high-altitude  nuclear  propagation  effects  are 
discussed.  The  results  are  timely  because  it  appears  that  there  is  a 
systematic  difference  between  (1)  the  power- law  index  that  characterizes 
the  phase  spectrum  of  radio  waves  that  have  propagated  through  a  striated 
barium  environment,  and  (2)  the  corresponding  power- law  index  for  most 
radio  waves  that  have  propagated  through  naturally  occurring  striations 
in  both  ionized  and  nonionlzed  media. 

It  is  generally  assumed  that  the  one- dimensional  in- situ  spectral 

density  function  (SDF)  that  characterizes  both  naturally  occurring  and 

-2±?  * 

barium  striations  has  the  form  K  ’  (Ossakow,  1979)*  For  neutral  tur- 

X  -5/3 

bulence  in  the  inertial  subrange,  the  Kolmogorov  theory  predicts  , 

whereas  the  nonlinear  steepening  of  Rayleigh  Taylor  and  E  X  B  coherent 

.2 

Fourier  modes  predicts*  .  The  number  of  published  in- situ  data  sets  is 

X  -2.3 

sparse;  however,  recently  analyzed  STRESS  data  suggest  *  for  barium 

striation. 

If  the  local  in-situ  characterization  of  the  irregularities  applies 
throughout  the  propagation  region,  then  the  corresponding  one- dimensional 
integrated- phase  SDF  has  an  index  one  larger  than  the  one- dimensional  in- 

-3+? 

situ  index.  Thus,  the  corresponding  phase  SDF  should  have  the  form  f  ”  * 
The  Wideband  satellite  data  have  consistently  shown  comparatively  shallow 
phase  SDFs  (f~p  with  2,2  <  p  ^  2.8)  at  all  latitudes  under  highly  varied 
propagation  conditions.  Similar  results  are  now  being  reported  for  inter¬ 
planetary  scintillation  and  line-of-sight  microwave  transmissions. 

If  one  accepts  the  simple  mapping  from  in-situ  to  phase  structure, 
there  is  a  discrepancy  between  the  in-situ  data  and  the  scintillation 
data.  Indeed,  the  scintillation  data  suggest  an  in-situ  index  smaller 

* 

References  are  listed  at  the  end  of  this  report. 
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than  the  Kolmogorov  value.  There  are  good  reasons  why  the  simple  mapping 
from  in-situ  to  phase  structure  might  not  be  invalid,  but  there  is  no 
question  that  it  is  the  effective  slope  of  the  integrated  phase  SDF  that 
determines  the  signal  structure. 

The  new  theoretical  results  we  have  alluded  to  show  that  there  is  a 
significant  change  in  the  signal  statistics  when  the  phase  spectral  slope 
approaches  or  exceeds  3.  Indeed,  the  propagation  modeling  task  is 
potentially  much  more  difficult  for  steeply  sloped  spectra.  In  light  of 
these  new  developments,  it  is  important  that  the  question  of  the  appro¬ 
priate  effective  spectral  index  for  propagation  modeling  be  resolved. 

To  introduce  the  problem,  consider  the  common  assumption  that  under 
conditions  of  sufficiently  strong  scatter,  the  complex  signal  v(t),  is 
Rayleigh  distributed.  This  means  formally  that  the  quadrature  components 
of  v(t)  are  independent,  identically  distributed,  gaussian  random  processes 
Since  the  second-order  statistics  completely  specify  a  gaussian  random 
process,  one  need  determine  only  the  mutual  coherence  function 

Rv(t,t')  =  <v(t)v*(t')>  (1) 

to  completely  characterize  the  temporal  structure  of  v(t)  at  a  single 
frequency.  In  particular,  all  higher- order  signal  moments  can  be 
computed  in  terms  of  Rv(t,t').  For  example, 

<I(t)I(t')>  -  <I(t)><I(t')>  =  lRv(t,t/)|2  (2) 

where 


effects  per  se  are  concerned,  departures  from  strict  Rayleigh  fading  are 
not  too  serious  (Johnson  and  Rino,  1979).  The  precise  value  of  the 
spectral  index,  however,  does  impact  predictions  of  time  structure,  fre¬ 
quency  coherence,  and  so  forth*  A  careful  assessment  of  the  signal 
structure  is  important,  therefore,  because  it  can  verify  the  value  that 
is  assigned  to  the  spectral  index* 

Indeed,  Rino  (1979b)  showed  that  in  a  three-dimensional  medium  with 
an  irregularity  spectral  density  function  (SDF)  of  the  form  Cgq  f 

converges  to  unity  in  the  strong- scat ter  limit  if  and  only  if  v  <  1*5* 
When  v  ^  1.5,  achieves  a  limiting  value  greater  than  unity,  which  can 
be  determined  from  the  formula 


2  _  4^2v  -  2 
4  5  -  2v 


(4) 


The  value  v  »  1.5  corresponds  to  a  one- dimensional  in-situ  SDF,  which, 

as  we  have  noted,  is  believed  to  be  associated  with  steepened  irregularity 
structures.  If  the  v  »  1.5  structure  model  is  assumed  to  have  global 
applicability,  then  Eq.  (4)  predicts  the  strong-scatter  limiting  value 
=  1.4.  Since  this  exceeds  unity  by  40%,  it  should  be  easily  detected 
in  any  statistically  valid  data  set.  Strong  focusing,  which  occurs 
before  the  limiting  value  is  achieved,  moreover,  will  cause  even  larger 
S^  values. 

In  addition  to  the  purely  phenomenological  impact  of  the  power- law 
slope,  there  are  ramifications  that  should  effect  predictive  codes. 

Indeed,  the  simple  form  that  is  being  used  for  the  phase  structure 
function  breaks  down  when  v  ^  1.5. 

What  is  happening  can  be  described  intuitively  as  follows:  As  the 
irregularity  SDF  steepens  (v  increases),  the  importance  of  the  omni¬ 
present  large-scale  structures  steadily  increases.  If  v  <  1.5,  these 
structures  manifest  themselves  only  in  the  first-order  moments  (e.g.,  as 
trend-like  phase  variations).  If  v  ^  1.5,  the  second-order  moments  (e.g., 
the  mutual  coherence  function)  develop  trend- like  components.  Finally, 
if  v  £  2.5,  even  the  fourth-order  moments  develop  trend- like  variations. 
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The  Wideband  satellite  data  clearly  favor  v  values  less  than  1.5  as 
evidently  do  interplanetary  scintillation  data,  and  data  from  light  and 
microwave  propagation  through  turbulent  atmospheres.  To  accommodate 
varying  v  values,  however,  a  change  in  the  usual  approach  to  propagation 
modeling  is  required.  Since  the  mutual  coherence  function  plays  a  central 
role  in  the  current  approach  to  predictive  modeling,  we  have  developed 
the  theory  in  detail  in  this  report.  The  development  shows  clearly  how 
inhomogeneities  manifest  themselves  and  what  their  impact  on  predictive 
modeling  is. 


II  THE  MUTUAL  COHERENCE  FUNCTION 


The  differential  equation  that  governs  the  development  of  the  mutual 
coherence  function  is 

7*4  — * 

dR  \P  *i  )P<i  '40  9  I 

dz  =  Zk  ^1  ~  "  *2^^15,:)2^  ^v^1,^2^  * 

2  2  2  2  2 

In  Eq.  (5),  =  d  /bpx^  +  S  /dpy^,  k  =  2jt/\,  and  is  the 

incremental  value  of  the  phase  structure  function, 

D(P1?P2)  =  (  [60(px )  -  6 0 ( P 2 )  ] >  .  (6) 

As  discussed  by  Rino  (1978),  all  the  moment  equations  derived  from 

the  parabolic  wave  equation  by  using  the  so-called  Markov  approximation 

contain  two  sets  of  terms.  One  set  of  terms,  which  involves  the  Laplacian 
2 

7^,  accounts  for  diffraction  effects  that  act  to  generate  small-scale 
structure  in  v(pf),  particularly  in  its  amplitude.  The  remaining  terms 
account  for  the  interaction  of  the  wavefield  with  the  randomly  irregular 
medium. 

Equation  (5)  is  written  in  its  isotropic  form.  The  general  form  is 
developed  in  Rino  (1978),  The  three-dimensional  structure  of  the  wave- 
field  in  an  anisotropic  medium  can  be  derived  by  straightforward  mani¬ 
pulations  of  the  isotropic  results  (Rino  and  Fremouw,  1977;  Rino,  1978, 
1979a).  Similarly,  the  temporal  structure  of  v(i?)  is  readily  derived  by 
using  an  appropriate  velocity  factor  to  convert  spatial  variations  to 
temporal  variations. 

As  a  general  rule,  whenever  the  diffraction  terms  impact  the  solution, 
the  particular  moment  will  depend  on  the  Fresnel  factor  Xz/4rt.  From 
Eq.  (5),  however,  it  is  readily  seen  that  if  D(p^,p2)  depends  on  the 
difference  variable  Ap  =  p 2  -  Pp  then  Rv(PpP2)  depends  only  on  Zp,  and 
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the  diffraction  terms  cancel*  It  follows  that  the  behavior  of  Rv(ppP2> 
depends  critically  on  the  spatial  homogeneity  of  the  structure  function* 

In  most  radio-wave  modeling  work  it  is  assumed  without  question  that 
any  departures  of  D(p^,p2)  from  strict  homogeneity  can  be  ignored.  To 
illustrate  that  there  is  good  reason  to  question  this  assumption,  consider 
the  signal  mean  (v(t))*  The  differential  equation  for  (v(t))  is 

V2<v(p)>  -  |r2\2  <£Ne(p»<v<£)>  .  (7) 

If  (£Ne(p))  is  strictly  independent  of  p,  then  (v(p))  is  constant.  In 
fact,  one  can  easily  obtain  the  solution  to  Eq.  (7)--namely, 

<v(p)>  =  exp  |~|<602>|  (8) 

where 


<S02>  =  r2X2 


G  l 


£<£n2>  d a 


(9) 


In  Eq.  (9)  G  is  a  purely  geometrical  factor  and  i  is  the  first  moment  of 
the  phase  spectral  density  function. 

Whether  or  not  there  is  a  spatial  regime  where  can  be 

considered  to  be  constant  can  only  be  determined  by  careful  data  analysis. 

Detrending  data  establishes  an  arbitrary  partitioning  between  large-  and 

small-scale  structures  such  that  the  detrended  component  exhibits  all  the 

2 

properties  of  a  homogeneous  process.  The  absolute  values  of  ($</>  )  and 
<v(p)>,  however,  are  determined  by  the  detrend  cutoff.  In  fact,  analysis 
of  the  Wideband  satellite  data  has  shown  that  even  with  detrended  data, 
(v(p))  is  not  well  behaved.  The  estimates  show  a  higher  degree  of  dis¬ 
persion  than  the  corresponding  estimates.  This  suggests  that  the 
diffraction  term  in  Eq.  (7)  cannot  be  ignored. 

Going  back  to  the  mutual  coherence  function,  if  D(p^,p2)  depends 
only  on  then  one  can  also  easily  obtain  the  solution  to  Eq.  (5)-- 
namely , 

8 


•-  •<  -jV  '  >*i.i**i*v 


(10) 


V&) =  exp  r 

where 


V*>  =  re^2 


i 


D(^p)  djj 


(ID 


We  note  that  the  initial  condition  =  1  at  z  =  0  is  implied  by  Eq.  (10). 
This  effectively  normalizes  the  average  signal  intensity  to  unity. 

2 

From  our  experience  with  the  first-order  moments  (&(f>  )  and  (v),  we 
take  it  as  a  guiding  principle  that  whenever  a  computed  signal  moment  of 
any  order  depends  strongly  on  qQ,  estimators  of  that  quantity  will  be 
dependent  on  how  they  are  measured.  The  effects  are  not  simply  due  to 
uncertainties  in  the  actual  value  of  qQ,  but  rather  to  diffraction  effects 
and/or  inhomogeneities  that  can  cause  significant  departures  from  the 
mathematical  limiting  form  of  the  particular  moment. 

To  pursue  this  further,  in  Section  III  we  discuss  the  mathematical 
model  that  is  being  used  for  the  phase  structure  function. 


9 


MhWW*  -lfc  r|M"  ■-  ■ 


Ill  THE  PHASE  STRUCTURE  FUNCTION 


To  simplify  the  analysis  without  loss  of  generality,  let  us  assume 
that  DG^,^)  does  not  vary  with  z.  For  a  path  of  length  Eq.  (11) 
then  becomes 

V(Pl’P2>  *  reX\  <t^e(Pl)  “  *  <12> 

The  integral  can  always  be  reintroduced  to  accommodate  slow  variations 
along  the  propagation  path.  The  structure  function  atuomatically  enters 
the  theory  as  discussed  in  Section  II.  It  has  been  recognized  for  some 
time,  however,  that  the  structure  function  intrinsically  suppresss  the 
low-frequency  content  of  the  field  it  is  performed  on. 

To  illustrate  this  directly,  consider  the  one-dimensional  form  of 
Eq.  (12).  It  is  readily  shown  that 

(  [A(x)  -  A(x  +  Ax)]2)  =  j"  [1  -  cos  (ouAx)]  $A(uj)  (13) 

where  $  (w)  is  the  power  spectrum  of  A(x).  Since  the  cosine  term  in 
a  2 

Eq .  (13)  varies  as  U)  for  small  u),  $^(uj)  ac  u>  ^  converges  independent 
of  any  explicit  outer-scale  cutoff  as  long  as  p  <  3.  This  is  exactly 
the  behavior  of  the  one-dimensional  phase  structure  functions  for  which 
p  =  2v  (Rino,  1979b).  It  might  also  be  noted  that  Eq.  (13)  is  identical 
in  form  to  an  MTI  (Moving  Target  Identification)  filter  that  is  used 
for  radar  clutter  suppression. 

A  convenient  SDF  model  that  is  both  mathematically  tractable  and 
accommodates  inner-scale  and  outer-scale  cutoffs  is 

$AN  (q)  =  CsqoV+1^T^F(^)  Kv+i(2sv/l  +  (q/qo)s)2e‘(V+^)/r(v-4)  (14) 

e  ~ 
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4.  ,■  >  ,*v ;  4  ■;  . 


where  e  =  q /q^  «  1  is  the  ratio  of  the  inner-  to  outer- scale  cutoffs. 

By  using  the  small- argument  approximation, 

Kviji(2x>  ~  3  r^v  +  ^  >ur'v  +  ^ .  (is) 

it  is  readily  shown  that  for  q  «  q^  and  e  «  1, 


$ 


AN 


e 


(q) 


c 

s 


(16) 


The  phase  autocorrelation  function  corresponding  to  Eq.  (14)  is 
R(y)  =  o^2v/es  +  (q0y/2)sV_^  kv_^(Vg3  +  (qcy/2)2)/N  (17) 


where 


N  =  2eV“*  K  a(2e) 

v-f' 

~  r(v  -  i)  . 


(18) 


The  approximation  in  Eq.  (18)  is  good  to  within  a  few  percent  for  v  <  0.5 
and  e  <  0.01.  In  the  small  e  approximation, 
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■  relV*  Wt 


r.,(y .-£> 

*r(\>4) 


-2v+l 


(19) 


It  is  notationally  convenient  to  let 


C 

P 


2,2, 
r  X  l 
e 


C 

P  s 


(20) 


so  that  the  phase  SDF  has  the  form  $  (q)  ~  C  q~(2v+l) 

2  P 

Section  II,  the  fact  that  depends  critically  on  qQ 
its  nonstationary  behavior. 


.  As  noted  in 
is  indicative  of 


11 


fTF 


The  form  of  the  structure  function  that  corresponds  to  Eq.  (17)  is 


Vy)  =  V(y)  =  CP 


2nT(y-i 


where 


£(y)  =  l 


+  (qoy/2)3' 


/ea  +  (qQy/2)3) /N 


In  Figure  1,  £(y)  is  plotted  versus  qQy  for  v  =  1.2  and  the  e  values 
0.01,  0.001,  and  0.0001.  It  is  readily  seen  that  the  impact  of  e  is 
negligible  even  for  comparatively  shallow  SDFs  where  its  effect  is  most 
pronounced.  For  v  =  1.5,  the  corresponding  curves,  which  are  shown  in 
Figure  2,  cannot  be  distinguished. 

Figure  3  shows  the  effect  of  varying  v.  At  first  glance  the 
variations  are  not  dramatic.  It  must  be  kept  in  mind,  however,  that 
the  significant  portion  of  the  curve  lies  to  the  left  of  yqQ  =  (2it)  *  = 
0.159,  since  this  corresponds  to  y  tf  ^  =  2it/q  .  As  v  increases,  the 
curves  are  displaced  to  the  right  of  y  =  indicating  that  structure 
is  being  determined  more  by  q^  and  less  by  the  power- law  slope. 

If  one  accepts  the  mathematical  form  of  Eq.  (22)  as  being  rigorously 

applicable,  the  £(y)  curves  admit  a  simple  interpretation.  Indeed,  from 

- 1  2 

Eq.  (10),  R  (4,)  =  e  when  crJ£(£,)  =  1.  It  follows  that  plots  of  £(y) 
v  a  p  d 

versus  yq  are  also  plots  of  l/a^  versus  l  ,q  ,  provided  that  <y  >  1. 
o  poo  p 

Since  £(f^)  saturates  at  unity,  it  is  obvious  that  no  meaningful  definition 

of  £,  can  be  obtained  for  small  <r  values.  The  fact  that  the  definition 
d  p 

of  depends  on  qQ,  however,  creates  problems  of  its  own. 

As  long  as  v  <  1.5,  there  is  a  simple  and  effective  way  to  deal  with 
the  latter  situation  that  has  been  used  in  turbulence  studies  for  over 
30  years.  It  is  shown  in  Appendix  B  that  for  sufficiently  small  qQ, 


Vy> 


FIGURE  3  PLOT  OF  NORMALIZED  STRUCTURE  FUNCTION  SHOWING  DEPENDENCE 
ON  SPECTRAL  SLOPE 


where 


9  c 
C2  =  -J2 
2  * 


2(1.5  -  v) 
r(v+0.5)(2v-l)2v_1 


is  the  phase  structure  constant.  To  demonstrate  the  validity  of  this 
asymptotic  approximation,  the  approximate  form  given  by  Eq.  (23)  is 
superimposed  on  the  exact  calculation  of  £(y)  in  Figure  4, 

The  plot  is  presented  on  an  expanded  scale  to  emphasize  the  range  of 
y  values  that  contribute  to  the  power- law  continuum- -viz,  y  <  i  [yqo  < 
0.159].  For  v  <  1,2,  the  small  approximation  works  essentially  over 
the  entire  y  <  &Q  range.  As  the  spectrum  steepens  (1.5  >  v  £  1.4),  the 
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FIGURE  4  PLOT  OF  SMALL  qQ  APPROXIMATION  TO  PHASE  STRUCTION  FUNCTION 
(valid  for  v  <  1.5)  SUPERIMPOSED  ON  EXACT  CURVES 


y  «  condition  becomes  more  stringent.  Thus,  as  long  as  v  <  1.4,  the 
structure-constant  approach  is  a  viable  one  for  the  entire  range  of  scale 
sizes  that  encompass  the  power- law  continuum. 

The  ultimate  test  of  any  model,  however,  is  how  well  it  reproduces 
data.  The  Wideband  satellite  data  have  consistently  shown  p  values  less 
than  3.  The  range  of  v  for  Wideband  is  1.2  ^  v  ^  1.3.  The  Kolmogorov 
value  is  v  =  4/3  =  1.333.  Excellent  results  have  been  obtained  in 

relating  the  intensity  coherence  time  under  strong-scatter  conditions 

2 

to  estimates  of  Ce  .  from  the  phase  data  (Rino  and  Owen,  1979).  In  the 
00 

analysis,  Eq.  (2)  was  effectively  used  to  relate  the  intensity  structure 
to  Equation  (2)  was,  however,  derived  by  calculating  the  strong- 

scatter  limiting  form  of  the  fourth-order  signal  moment. 


Attempts  to  measure  the  mutual  coherence  function  directly  have 
shown  that  Eq.  (10)  is  strictly  valid  for  a  limited  range  of  y  values. 

For  example,  contrary  to  the  predictions  of  Eq.  (10),  measured  values  of 
(v(t)v*(t/))  invariably  go  negative.  At  the  present  time,  attempts  are 
being  made  to  average  several  estimates  together,  but  it  appears  that 
the  departures  from  the  idealized  model  are  not  simple  random  fluctu¬ 
ations  that  ultimately  average  out  to  the  "correct"  value. 

What  remains  is  to  consider  the  situation  for  v  £  1.5.  As  noted  in 
Section  I,  the  radio- wave  propagation  community  has  accepted  the  v  =  1.5 
value  as  a  norm,  whereas  the  Wideband  satellite  data  have  consistently 
shown  smaller  values.  The  only  carefully  analyzed  data  that  have 
indicated  steeply  sloped  SDFs  have  come  from  in-situ  equatorial  rocket 
probes  and  in-situ  as  well  as  propagation  data  from  the  STRESS  barium 
exercises.  The  early  equatorial  rocket  probes  can  be  dismissed  since 
they  only  penetrated  bottomside  spread  F.  The  preliminary  results  from 
the  recent  successful  Kwajalein  rocket  campaign  show  that  bottomside 
spread-F  is  a  mere  ripple  when  compared  with  well  developed  topside 
structures . 

A  barium  cloud  is  a  highly  localized  structure  whose  evolution 
can  be  optically  tracked.  A  detailed  analysis  of  the  STRESS  barium 
cloud  ESTER  by  MRC  has  shown,  firstly,  that  structures  larger  than  1  km 

cannot  be  reliably  included  in  the  power- law  continuum,  and  secondly, 

-2 

that  the  smaller-scale  structures  fall  off  more  steeply  than  K 

The  only  radio  propagation  data  were  obtained  from  UHF  transmissions 
that  were  highly  disturbed.  ESL  has  used  a  ,rback  propagation"  technique 
to  reconstruct  the  integrated  phase  SDF,  and  the  preliminary  results 
seem  to  be  in  agreement  with  the  in-situ  data  in  that  the  phase  SDF  falls 

-3 

off  more  steeply  than  f  •  Under  such  conditions,  however,  one  should 
expect  significant  departures  from  Rayleigh  statistics.  The  small  size 
of  the  largest  structure  may  make  tests  of  the  statistics  difficult.  A 
straightforward  test  of  the  back- propagation  procedure,  however,  is  to 
measure  the  mutual  coherence  function  directly  and  see  how  it  conforms 
to  model  calculations  based  on  i* (y )  for  v  ^  1.5. 


One  would  like  to  believe  that  the  barium  data  accurately  reproduce 
what  it  happening  locally.  Over  much  larger  regions,  however,  the 
integration  over  many  structured  subregions  may  well  act  to  produce  a 
shallower  integrated  phase  spectral  slope.  In  that  case,  the  effective 
v  index  for  calculating  propagation  effects  is  less  than  1.5  and  the 
diffraction  theory  is  greatly  simplified. 

Recent  analysis  of  30-GHz  scintillation  data  over  long  atmospheric 
paths  has  shown  phase  spectral  slopes  significantly  lower  than  the 
expected  8/3  value  based  on  the  Kolmogorov  theory.  A  random  patch  model 
has  been  invoked  to  reconcile  the  discrepancy.  The  model  is  supported 
by  high-resolution  radar  data  showing  that  clear-air  turbulence  develops 
in  distinct  narrow  regions. 

However  the  issue  is  resolved,  there  is  a  significant  discrepancy 
between  structure  models  inferred  from  the  Wideband  satellite  data  and 
the  structure  models  that  are  being  extracted  from  the  STRESS  barium 
data.  The  recent  Kwajalein  rocket  data  will  undoubtedly  shed  some  light 
on  the  issue,  but  it  must  be  resolved  if  a  viable  predictive  propagation 
code  is  to  be  developed. 
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IV  THE  TAYLOR  SERIES  APPROXIMATION 


Several  authors  have  applied  a  Taylor  series  approximation  to  £(y) 
to  simplify  analyses  involving  integral  expressions  containing  various 
combinations  of  structure  functions  (e.g.,  Taylor  and  Infosino,  1976; 
Buckley,  1971).  As  discussed  by  Rino  (1979b),  however,  these  series 
converge  very  slowly,  particularly  in  the  small  limit.  In  a  recent 
analysis  of  sound  propagation  in  sea  water,  Dashen  (1977)  concluded  that 
the  Taylor  series  approximation  cannot  be  used  in  power- law  medium. 

Dashen,  however,  was  concerned  with  SDF  corresponding  to  v  <  1.5, 
and  the  Taylor  series  approximation  improves  as  v  increases.  Indeed, 
this  is  expected  because  the  derivatives  involved  are  equivalent  to 
moments  in  the  spectral  domain.  The  moments  are  increasingly  large 
unless  the  SDF  falls  off  very  rapidly.  An  inner-scale  cutoff  keeps  the 
derivatives  finite,  but  the  convergence  of  the  series  is  still  very  slow 
if,  indeed,  proper  convergence  even  occurs. 

In  any  case,  the  formal  Taylor  series  expansion  of  £(y)  is  easily 
derived.  If  we  use  the  notation 

Ri(z)  =  2zV"^Kv^(2z)/N  (25) 

and 


f(y)  =  Jp  +  (yqQ/2)3 


(26) 


it  is  shown  in  Appendix  A  that 


*C(y> 

3y 


-  <Pv.3/2  (£(y»  v« 


(27) 


From  Eq.  (27)  all  higher-order  derivatives  can  be  generated.  For  example 


(28) 


’  -*v-5/2  ’o^4  *  *v-3/2  <f  1 <*»  1o/2  • 

by 


In  the  y  =  0  limit,  f(y)  =  e,  and  only  the  even  derivatives  survive. 
The  formal  Taylor  series  can  be  written  as 


i(y)  =2  Dn(qoy)2n  (29) 

n=l 


where 


D 

n 


a2.n  £h±. 

(2n)!&y2n 


(30) 


The  coefficients  are  derived  in  Appendix  A.  For  our  purposes  we  need 
consider  only  the  first  coefficient, 


D 


_  1  v-3/2 
1  2e 


Kv-3/2(2e)/N 


(31) 


n3/2-v)  3 

4r(v4) 


log(2e) 

2r(v4) 


v  <  3/2 


(32) 

v  —  3/2  . 


From  Eq.  (32)  we  see  that  for  v  <  1.5,  is  very  sensitive  to  e. 

In  fact,  the  higher-order  derivatives  are  even  more  sensitive  to  e  and 
the  Taylor  series  approach  is  useless,  as  noted  by  Dashen  (1977).  It 
appears  that  the  series  is  divergent  except  for  very  small  y  values. 

The  dependence  on  e  decreases  as  v  increases.  If  v  >  1.5,  does  not 
depend  on  e,  but  higher  derivatives  do,  and  convergence  problems  persist. 

In  Figure  5  the  quadratic  approximation  is  shown  for  v  =  1.2,  1.5, 
and  2.0,  The  approximation  Is  acceptable  for  v  =  1.5,  provided  that 
y  <  It  improved  as  v  exceeds  1.5,  but  the  series  approximation 

cannot  be  used,  as  we  noted  previously.  We  emphasize  that  neither  the 
asymptotic  nor  the  quadratic  approximation  admits  y  values  larger  than 
lQ  =  2rt/qo. 
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FIGURE  5  PLOT  OF  QUADRATIC  APPROXIMATION  TO  PHASE  STRUCTURE  FUNCTION 
SUPERIMPOSED  ON  EXACT  CURVES 


V  SUMMARY  AND  DISCUSSION 


In  this  report  we  have  summarized  some  new  results  in  propagation 
theory  (Rino,  1979a, b,c;  Rino  and  Owen,  1979),  and  have  discussed  their 
ramifications  for  predictive  modeling.  The  theory  was  motivated  by  the 
Wideband  satellite  data  which  clearly  showed  the  inhomogeneous  structure 
of  the  integrated  electron  density.  One  observes  slow  trend- like  vari¬ 
ations  in  phase  because  of  the  dominance  of  large-scale  structures  in  a 
power- law  environment. 

The  large-scale  structures  potentially  impact  all  signal  moments. 

The  impact  can  be  nil  or  totally  dominating,  depending  on  the  order  of 
the  moment  and  the  power- law  index.  The  first-order  moment  of  the 
complex  signal  is  exponentially  dependent  on  the  integrated  phase;  thus, 
it  exhibits  the  same  nonstationary  behavior  as  the  integrated  phase. 

The  second- order  complex  signal  moment  (the  mutual  coherence 
function),  which  is  used  directly  or  indirectly  for  predicting  essentially 
all  strong-scatter  propagation  effects,  is  not  sensitive  to  the  large-scale 
structures  if  the  effective  spectral  index  v  is  less  than  1.5,  which  is 
evidently  the  case  for  all  naturally  occurring  irregularities.  If 
v  >  1,5  as  the  STRESS  data  seem  to  indicate,  the  large-scale  structures 
dominate  the  second-order  moments.  Predictive  modeling  then  depends 
critically  on  the  outer-scale  parameter. 

The  Rayleigh  hypothesis  breaks  down  when  v  >  1.5.  Insofar  as 
system  effects  per  se  are  concerned,  departures  from  the  gauss ian  model 
that  is  invariably  used  for  analyzing  or  simulating  disturbances  are  not 
too  important.  For  predictive  modeling  of  these  parameters  that  charac¬ 
terize  the  signal  structure,  however,  there  is  a  significant  impact. 

Because  of  this  fact  we  have  suggested  some  independent  checks  of  the 
back- propagation  analysis  of  the  STRESS  scintillation  data.  Specifically, 
the  departures  from  strict  Rayleigh  statistics,  in  particular  >  1, 
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that  must  accompany  steeply  sloped  SDFs  should  be  verified,  if  possible. 
The  mutual  coherence  function  itself  should  be  carefully  measured  to 
verify  that  it  conforms  to  model  predictions. 
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Appendix  A 


THE  TAYLOR  SERIES  APPROXIMATION 
TO  THE  PHASE  STRUCTURE  FUNCTION 


The  normalized  phase  structure  function  can  be  written  as 

■r(y)  =  i  -  ^(f(y»  (a-  1 ) 


where 


^v-i(z)  =  2zV"^ 

kv4(2z)/N 

(A- 2) 

N  =  2gV^ 

Kv-i(2e) 

(A- 3) 

and 


f(y)  =  Ve5  +  (qQy/2)s 


(A-4) 


It  is  clear  from  Eq.  (A-4)  that  is  simply  a  scale  factor.  We  shall 
derive  a  formal  Taylor  series  representation  for  Eq.  (A-l)  with  the  aim 
of  obtaining  a  simple  approximation  to  £(y)  that  is  valid  for  e  «  1,  and 
as  large  a  yq^  range  as  possible. 

To  start  with,  consider  the  first  derivative, 


itM 

dy  by 


(A-5) 


The  derivatives  of  products  such  as 
easily  shown  that 


K^^(2z)  are  known, 


and  it  is 


f(y) 


-  (*<y»  v/2 


(A- 6) 


A-l 


Using  Eq.  (A- 6)  we  have 


^  ■  *v-3/2  1»y/2  <A-7) 

’Z?"1  ■  ■pJ-5/2|f(y)|  qoy2/4  +  3/2  ,f(y)*  <<o/2  <A'®> 

5y 

and 

3 

^3^  =  ^v- 7/2  (f(y))  V3/8  ‘  ^v-5/2  <f(y)>  qo3y/4  *  (A_9) 

Sy 

The  n^  derivative  of  £(y)  has  Nq  =  (n-l)/2+l  terms  if  n  is  odd,  and 
N0  =  n/2+1  terms  if  n  is  even.  The  first  term  always  has  the  form 

(A-10) 

where 

c"  =  -  1/2  c"-1  .  (A- 11) 

We  let  denote  the  coefficient  of  the  i^  term  in  the  n^  derivative. 

If  n  is  odd,  the  terms  i  =  2,  3,  . . . ,  Nq  have  the  form 
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\' 
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^(v-l/2)-n+(i-l) 


2(n-i+l)  (n—  2  (  i—  1 ))  n 
qo  y  ^i 


where 


c[  =  [(n+3-2i)  |cj;j| 


(A- 12) 


(A- 13) 


The  same  formula,  Eq.  (A- 12),  applies  to  the  even  derivatives,  with  the 
exception  that  the  last  coefficient,  / 2-hl 5  et^ua^  t0  t^ie  previously 
computed  last  odd  coef f icient--i.e. , 
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(n  even) 


(A- 14) 
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n/2+1 


«  C 


n- 1 
n/2 


The  formal  Taylor  series  can  be  written  as 


£(y)  Dn(yqo)2n  (A_15) 

n=l 


where 


D 

n 


a2n 

(2n)J 


£1x1 


y~o 


(A- 16) 


th 

Only  the  Ng  even  term  gives  a  finite  contribution.  Hence,  from  Eq. 
(A- 12) 


Dn  *  (2n) I  /?(v-l/2)-n  f(y)  ly=0  Cnfl 

1  v-l/2-n  /o  \  /m  /a  i  7 \ 

(2n) !  e  KV-l/2-n(2e)/N  Cn+1  *  (A-17) 

Note  that  as  long  as  v-l/2-n  >  0,  Dn  is  nearly  independent  of  e  for 
small  c •  For  large  n  values,  oc  e^v  ^  ^n,  which  becomes  arbitrarily 
large  with  increasing  n.  It  follows  that  the  range  of  applicability  of 
the  Taylor  series  approximation  is  very  restrictive.  Recent  analysis  by 
Dashen  (1977)  suggests  that  it  can  only  be  used  for  propagation  environ¬ 
ments  with  a  sharp  cutoff. 
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Appendix  B 


THE  SMALL  q  APPROXIMATION 
o 

Since  (j.  «  q"2v+1,  it  follows  that  <j2  (1  -  JC(y ))  converges  to  the 
©  o  © 

indeterminate  form  —  when  yq^  «  1  and  q^  is  small.  If  v  <  1.5,  however, 
o  o  o 

L’Hopital's  rule  can  be  applied  to  compute  the  limiting  form  of 

D(y)  =  cr?  (l  •  £(y))«  We  first  let  e  -*  0  to  simplify  £(y)  to 
0 

<E(y)  ~  1  -  2|qoy/2|V_1/2Kv_1/2(qoy)/r(v-l/2)  (B-l) 


which  is  the  form  used  in  Rino  (1979b).  The  quantity  of  interest  is 


lim  £(v) 

q  — *0  2v-l 
° 


(B-2) 


Using  the  results  of  Appendix  A  to  calculate  the  derivative  of  £( y), 
from  L'Hopital's  rule  we  have 


lim  JE(y) 
q  — *0  2v- 1 

o 


lim  q0(v/2)V"3/2  Vv^v^y2 


v° 


r(v-l/2)(2v-l)q 


2v-l 


(B-3) 


As  long  as  v  <  1.5,  the  small- argument  approximation  of  3/2^0^  shows 

that  the  q  terms  cancel  and  we  have 
o 


lim 
q  -*0 


a‘ *6  (y) 


C 

=  _E 

2* 


2(1. 5-v) 


r(v+0.5)(2v-l)2 


2v-l 


,y,2-1 


(B-4) 


This  result  has  been  used  in  studies  of  neutral  turbulence  for  decades • 
The  coefficient  that  multiplies  |y|^V  *  is  the  phase  structure  constant* 

In  studies  involving  radio-wave  interactions  with  turbulent  neutral 
atmospheres,  only  the  structure  function  enters  the  problem.  Indeed,  the 
structure  constant  is  typically  modeled  directly.  Only  recently  has 


B-l 


erious  consideration  been  given  to  propagation  ettects  in  environments 
haracterized  by  steeply  sloped  power- law  SDFs* 
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